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ABSTRACT 

Dust in galaxies forms and evolves by various processes, and these dust processes change the 
grain size distribution and amount of dust in the interstellar medium (ISM). We construct a 
dust evolution model taking into account the grain size distribution, and investigate what kind 
of dust processes determine the grain size distribution at each stage of galaxy evolution. In 
addition to the dust production by type II supemovae (SNe II) and asymptotic giant branch 
(AGB) stars, we consider three processes in the ISM: (i) dust destruction by SN shocks, (ii) 
metal accretion onto the surface of preexisting grains in the cold neutral medium (CNM) 
(called grain growth), and (iii) grain-grain collisions (shattering and coagulation) in the warm 
neutral medium (WNM) and CNM. We found that the grain size distribution in galaxies is 
controlled by stellar sources in the early stage of galaxy evolution, and that afterwards the 
main processes that govern the size distribution changes to those in the ISM, and this change 
occurs at earlier stage of galaxy evolution for a shorter star formation timescale (for star 
formation time-scales = 0.5, 5 and 50 Gyr, the change occurs about galactic age t ^ 0.6, 2 
and 5 Gyr, respectively). If we only take into account the processes which directly affect the 
total dust mass (dust production by SNe II and AGB stars, dust destruction by SN shocks, and 
grain growth), the grain size distribution is biased to large grains (a ^ 0.2-0.5 /im, where a 
is the grain radius). Therefore, shattering is crucial to produce small (a < 0.01 /xm) grains. 
Since shattering produces a large abundance of small grains (consequently, the surface-to- 
volume ratio of grains increases), it enhances the efficiency of grain growth, contributing to 
the significant increase of the total dust mass. Grain growth creates a large bump in the grain 
size distribution around a ^ 0.01 /im. Coagulation occurs effectively after the number of 
small grains is enhanced by shattering, and the grain size distribution is deformed to have a 
bump at a ~ 0.03-0.05 /im at < ^ 10 Gyr We conclude that the evolutions of the total dust 
mass and the grain size distribution in galaxies are closely related to each other, and the grain 
size distribution changes considerably through the galaxy evolution because the dominant dust 
processes which regulate the grain size distribution change. 

Key words: dust, extinction - galaxies: evolution - galaxies: ISM - ISM: clouds - galaxies: 
general - stars: formation 



1 INTRODUCTION 



Dust is one of the most important factors for the understanding of galaxy evolution. Since hydrogen molecules are efficiently formed on the 
surface of dust grains, the molecular formation rate is much larg er than the case without d ust. Such an enrichment of molecular abundance 
by dust realizes a favorable condition for star formation (e.g.. iHirashita & Ferrara 20021) . Dust grains also absorb stellar light mainly at 
ultra violet and opt ical wavelengths and re-emit in the infrared. Consequently, dust affects the spectral energy distribution (SED) of galaxies 
C^g.- IXakagi et alf2 003). Furthermore, the formation rate of hydrogen molecules on the grain surface and th e mass absorption coefficient of 
radiation depend strongly on the grain size distribution (e.g., Hirashita & Ferraral2002 : Takeuchi et al ]|2003h . 
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Dust grains form by condensation of elements heavier than helium (i.e., metals). Metals are mainly supplied from asymptotic giant 



branch stars (AGB stars) and supernovae (SNe), and part of them condense into dust gr ains (e.g,. 



Mathis 



suppl i ed by stars but are also destroyed by SN shocks in the interstellar medium (ISM) (e.g., Jones. Tielens & HoUenbachl 19961 : 



Dust grains are not only 



Nozawa et al 



20061 : IZhukovska. Gail. & Tr ieloflll2008l) Furthermore, it is thought that metal accretion onto the surface of grains in the ISM (referre d to 



as "gr ain growth" in this paper) is an important process for explaining the amount of dust in the Milky Way (e.g.. Draine 2009: Pipino et alj 
201lh . To the present day, the re have been a lot of studies that i nvestigate the ev ol ution of the total dust mass in galaxies by taking into 
account these processes (e.g ., 'Dwek & Scalo"l98o'; 'Dwek"l998'; 'Hiras hitj[l999al lbl: llnoudboosl boul : Icalura. Pipino & Matteuccil l2008l : 



dust destruction £ind grain growth depend on the grain size distribution. Thus, we should consider the evolution of the grain size distribution 
to understand the total dust mass precisely. 

The grain size distribu tion is derived from observed extin ction curves (which mainly depend on the grain size distribution and the 
grain species). According to lMathis. Rumpl & Nordsieckl ( Il977n . if spherical grains are assumed, the extinction curve in the Milky Way is 
reproduced by f{a)da oc a~"^'^da (0.005 < a < 0.25 ^m; this grain size distribution is referred to as the MRN distribution), where a is the 
grain radius and f{a)da is the number density of grains in size interval [a, a + da] (seel Kim, Martin. & Hendrvll 19941 : IWeingartner & Draind 
for more detailed fitting to the Milky Way extinction curve). The situation seems to be very different for distant galaxies. Recently, 



2001 



Gallerani et al. i201ol) discussed the extinction curves of seven quasars at high redshift (3.9 ^ z ^ 6.4). They showed that these extinction 



curves tend to be flat at wavelengths < 0.2 fim in the quasar's rest frame. The difference between extinction curves in distant and nearby 
objects may indicate that different processes dominate the dust evolution at different epochs. 

In y oung galaxies. Type I I SNe (SNe II) are thought to be the dominant sources of dust because of short lifetime of their progenitors. 
However, Valiante et al. 1 2009 1) suggested that AGB stars are also important sources of dust production even at galactic age less than 1 Gyr. 
In addition, g rain growth is expected to be th e dominant process to increase dust mass in galaxies if the metallicity becomes larger than a 
certain value i Inoue l201ll : lAsanoetal.ll2013l) . Furthermore, if the metallicity reaches a sub-solar value, gr ain-grain collisions in the ISM 
(shattering and coagulation) become efficient enough to change the grain size distribution significantly (e.g.,|Hirashita & Yanll2009h . We call 
all processes affecting the grain size distribution 'dust processes'. 

These dust processes affect the different sizes of grains in galaxies. iNozawa et alj ( 120070 showed that SNe II supply relatively large 
grains (a > 0.01 /xm) into the ISM because small grai ns are destroyed by reverse shocks before they are ejected into the ISM (see also 
Bianchi & Schneidenl2007l :ISilvia. Smith, & Shull 2012). The size distribution of grains produced by AGB stars is thought to be biased to 



large (^ 0.1 /xm) sizes (e.g., | Winters et al. 1 997. : , Groenewegen. 1997 ; Yasuda & Koz asa 2012). Furthermore, the smaller grains in the ISM 
are more easily destroyed by interstellar shocks driven by SNe (e.g., Nozawa et al. I2OO6I) . If grain growth occur s, since the timescale o f 
this process is proportional to the volume-to-surface ratio of a grain, smaller gr ains gro w more efficiently (e.g., Hirashita & Ku j 201 1 ). 
After the dust grains are released into the diffuse ISM, shattering can also occur. lYan et a l. (20_04) showed that large grains (a > 0.1 fim) 
acquire larger velocity dispersions than the shattering threshold velocity if the grains a re dynamically coupled with magnetized interstellar 
turbulence. Sh attering is indeed a promising mechanism of small-grain production (e.g., l Hirashital20 l"q). Shattering also occurs in SN shocks 
I Jones. Tielens & Holienbach 1996). In dense regions, coagulation can occur, so that the grain size distribution shifts to larger sizes (e.g., 
Ormel et al .boogtlffirashita & Yanl2009l) . The various dust processes above in galaxies occur on timescales dependent on the metallicity, the 



total dust amount, the grain size distribution, and so on. Hence, it is crucial to consider all dust processes in a unified framework to understand 
the evolution of both the total dust amount and the grain size distribution. 

There have been a number of studies on the evolution of the grain size distribution in galaxies. iLiffman & Clavton ( Il989h discussed the 
evolution of grain size distribution consi dering dust destruction by SN shocks and grain growth. However, they did not consider shattering 
and coagulation by grain-grain collisions. O'Dormell & Mathi3 ( 1997 ) suggested a dust evolution model in a multi-phase ISM [warm neutral 
medium (WNM) and cold neutral medium (CNM)], and also considered the coUisional p rocesses of dus t grain s. However, they did not 
consider the size distribution of grains released by stars in order to simplify their model. iHirashita et al. ( 20I0l) discussed the grain size 
distribution in young starb urst galaxies. They assum ed that SNe II are the source of dust in these galaxies and focused on the production of 



small grains by shattering. 



Yamasawa et al 



( 20Ilh constructed a dust evolution model taking into account dust formation and destruction 
by SNe II along with the formation and evolution of galaxies. However, since they focused on galaxies in the high-z Universe, they did not 
consider dust formation in AGB stars, grain growth, shattering and coagulation. 

In this work, we construct a dust evolution model taking into account the dust formation by SNe II and AGB stars, dust destruction by 
SN shocks, grain growth, and shattering and coagulation, to investigate what kind of dust processes determine the grain size distribution at 
each stage of galaxy evolution. In our model, we do not consider mass exchange among various ISM phases in detail (e.g., Ikeuchi & Tomital 
1983h . but our results contain the contributions of dust processes in the two ISM phases, WNM (~ 6000 K, 0.3 cm ) and CNM (~ 100 K, 



30 cm ) by assuming these mass fractions in the ISM to be constant. 

This paper is organized as follows: in Section |2] we introduce the dust evolution model based on chemical evolution of galaxies. In 
Section [3] we examine the contribution of each dust process to the grain size distribution. Section |4] is devoted to the discussion on what 
kind of dust processes regulate the grain size distribution in galaxies. We conclude this work in Section |5] Throughout this paper the solar 
metallicity is set to be Zq = 0.02 dAnders & Grevesselll989h . 
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2 GALAXY EVOLUTION MODEL 

In this Section, we introduce our dust evolution model in a galaxy. First, we show the basic equations of the chemical evolution model. We 
then describe the dust evolutions based on the chemical evolution model, involving dust production by SNe II and AGB stars, dust destruction 

by SN shocks, grain growth, and shattering and coagulation by grain-grain collisions. 

Some grain processing mechanisms work in a different way in a different ISM phase ( O'Donnell & MathisI 1997 ). In this work, while 



we use a one-zone model to examine the representative properties of a galaxy, we consider the effects of the dust processes in WNM and 
CNM by introducing the mass fractions of WNM and CNM, 77WNM and 77cnm- Considering temp eratures less than 10 '^ K, we find that an 
equilibrium state of two thermally stable phases (WNM and CNM) is established in the ISM (e.g.. lWolfireetalJl2003h . Thus, we calculate 



dust evolution taking into account a two-phase neutral ISM. We also assume that the galaxy is a closed-box; that is, the total baryon mass 
Mtot (the sum of the stellar mass and the ISM mass in the galaxy) is constant. Since i\ftot is just a scale factor in our work, the total dust 
mass just scales with A/tot - Throughout this paper Aftot is set to be 10^" Mq. 

Inflow and outflow are not considered in our model for simplicity. Since inflowing gas is considered to be not only metal poor but also 
dust poor, the abundance of both metals and dust are diluted with the same (or similar) fraction by inflow. This effect is degenerate with a 
slower chemical enrichment, under a longer rsp, where rsp is the star formation timescale. As for outflow, since ISM components (namely 
gas, metals, and dust) are blown out of a galaxy, the total gas mass in a galaxy decreases. In this case, star formation rate decreases at earlier 
phase of galaxy evolution; that is, the effect of outflow is degenerate with a shorter rgp. Thus, we just absorb the effects of inflow and outflow 
into rsp. 



2.1 Chemical evolution model 

In this subsection, we describe our model of chemical evolution in a galaxy. From the above assumptions, the equations of time evolution of 
the total stellar mass, Af*, the ISM mass. A/ism, and the mass of a metal species X, A/x, in the galaxy are expressed as 

dM,{t) 



dt 

dAfiSNift) 
dt 
dMxit) 



= SFR{t)-R{t), (1) 
= -SFR{t) + R{t), (2) 
= ^Zx{t)SFR{t) + Yx{t), (3) 



dt 

where t is the galaxy age, SFR(t) is the star formation rate, Zx = Afx/AfisM, and R{t) and Yx{t) are the masses of the total baryons and 
total metal species X released by stars in a unit time, respectively. In this paper, we consider two dust species, carbonaceous dust and silicate 
dust, and we adopt two key elements of dust species (X = C for carbonaceous dust and X = Si for silicate dust) in calculating dust evolution 
(see Section [Z2l for details). We adopt A/, (0) — 0, A/ism (0) = Aftot, and A/x(0) = as initial conditions. 

In our work, we adopt the Schmidt law for the SFR: SFR oc M^-^ jschmidtil959i) . and the index n is thought to be 1-2 observationally 
(e.g., Kennicutill 19981 ). We here adopt n = 1, 



SFR(t) = Ml^, (4) 

TSP 

where the star formation timescale tsf is a constant. For comparison, the case with n — 1.5 is also shown in Appendix I All As long as we 
adopt the same star formation timescale at f = 0, there is little difference between the two cases with n = 1 and 1.5. 
R{t) and Yx{t) are written as 



100 M 



O 



R{t) = I [m-u;{m,Z{t-T,„))](t>{m)SVR{t)dm, (5) 

it(t) 



(.100 Mq 

Yxit) = / mxiTn,Z{t-r„,))(t){m)SFR{t)dm, (6) 

Jmcut(t) 

where </>(m) is the stellar initial mass function, Tm is the lifetime of a star with mass m at the zero-age main sequence, Z is the metallicity 
(= T^xMx/Mism), and uu{m, Z) and mxijn, Z) represent the mass of remnant stars (white dwarfs, neutron stars or black holes) and the 
mass of metal species X ejected b y a star of mass m and metallicity Z, respectively. For the lifetime of stars, we adopt th e formula derived b y 



Raiteri. Villata & Navarro 



el ected b y a star ot mass m and metallicity Z, respectively, tor tne liietime 01 stars, we adopt tn e tormula derived b y 
3 (ll996l) . and the formula is obtained by the fitting to the stellar models of the Padova group jBertelli et alJI 19941) . 



Since its metallicity dependence is weak, we always adopt the stellar lifetime for solar metallicity as a representative value. The lower bound 
of the integration, mcut(t) is the mass of a star with Tm = t. We adopt the SalpeterlMF for stellar mass range 0.1-100 Mq (iSalpeter.19551) : 

(j>{m) oc mT'' , (7) 

where q is set to be 2.35, and the normalization is determined by 

(.100 Mq 

/ mcj>{m)dm = 1. (8) 

JO.l Mq 



© 2013 RAS, MNRAS 0Q0.[Tll20l 



4 Asano et al. 



Table 1. parameters for each dust species 



Species 


X 


5X 


mx [amu] 


s [g cm ^] ^ 


fshat [km s 1] 


7 [erg cm ^] 


£ [dyn cm-'-^] 




Graphite 


C 


1.0 


12 


2.26 


1.2 


75 


1.0 X 10" 


0.32 


Silicate 


Si 


0.166=* 


28.1 


3.3 


2.7 


25 


5.4 X 10" 


0.17 



Note. X is the key element of dust species, gx is the mass fraction of the key element X in the grains, mx is the atomic mass 
of X, s is the bulk density of dust grains, iighat is the shattering threshold velocity. 7 is the surface energy per unit area of 

grains, E is Young's modulus, and u is Poisson' s ratio. 

^We assume Mgj^ ^Fcq. 98104 for the composition of silicate iDraine & Leell984l) . 

'HDraiiie& Led 11984. 

HChokshi. Tielens. & HolleiibacHll993l l. 



To check the variation of the results with q, we examine the case with q = 1.35 (a top heavy IMF) in Appendix I A2I For q — 1.35, the 
processes in the ISM occur at earlier phases of galaxy evolution than for q = 2.35, because a larger amount of dust is supplied by stars. 
However, we find that the sequence of the dominant dust processes along the age does not change so the following discussions are not altered 
significantly by the change of q. Thus, we only consider q — 2.35 in the following discussion. 

To calculate Eqs. ^ and (|6}, we quote the remnant and metal mass data of stars with mass m and metallicity Z from some previous 
works. We assume that the mass ranges of AGB stars and SNe II are 1-8 M0 and 8-4 Mp), re spectively, and that all stars with initial masses 
more than 40 Mq evolve into black holes without ejecting any gas, m etals or dust 1 Heger et al. 2003) . The data for AGB stars with mass 
1-6 Mo and metallicity Z = (0.005, 0.2, 0.4, 1.0) Zb is taken fromkarakasl ilOld) and thedata for SNe II with mass 13-40 Mq and 



metallicity Z — (0.0, 0.05, 0.2, 1.0) Zq is from Kobavashi et alj (200a). We interpolate and extrapolate the data for all values of mass and 
metallicity (also for the dust data in Sections [2.2.1 I and l2.2.2b . 



2.2 Dust evolution 



For dust evolution, we consider dust production by SNe II and AGB stars, dust destruction by SN shocks in the ISM, grain growth in the 
CNM, and shattering and coagulation by grain-grain collisions in the WNM and CNM. In this work, as mentioned in Section[T] we assume 
a two-phase ISM (WNM and CNM) to calculate the dust evolution (see also Section [2.2.5b . 

We neglect the contribution of Type la SNe (SNe la) to the production of metals and dust, and the destruction of dust. iNozawa et alj 



1 201 ih showed that SNe la release little dust into the ISM. Furthermore, dust de struction by SNe la is expected to be insignificant to the 



total dust budget in galaxies (less than 1/10 of the contribution of SNe II; Calura. Pipino & Matteucci (^08)). As for metals, although 



Nomoto et al. 



1 19970 showed that the contribution of SNe la to the silicon and car bon enrichment in the ISM can be comparable to that of 



SNe II, the ratio between SN la rate and SN II rate is unknown | Nomoto et alj 1 1997 ) suggested that it is about 0. 1 taking into account a chem- 



ical evolution model]. Thus, to simplify the discussion, we neglect the contribution from SNe la, keeping in mind a possible underproduction 
of metallicity. 

The dust production data we adopt c ontain a lot of dust specie s (C, Si, Si02, SiC, Fe, FeS, AI2O3, MgO, MgSiOg, FeSiOa, Mg2Si04, and 
Fe2Si04) (Nozawa et al. 2007 : Zhukov ska. Gail. & Trieloffl20oi) . However, the physical properties of grain species other than carbonaceous 



and silicate grains are not fully known. Hence, we categorize all grain species other than carbonaceous dust as silicate and calculate their 
growth, shattering, and coagulation by adopting the physical parameters of silicate grains. In particular, after grain growth and coagulation 
occur, the dust species categorized as silicate dust do not evolve separately and our simplification can avoid the complexity arising from 
the compound species. In fact, the mass of dust grains ejected by SNe is dominated by Si grains, which would grow into silicate grains 
in the oxygen-rich environments such as molecular clouds. For carbonaceous dust, we adopt material properties of graphite. The adopted 
parameters of these two grain species are shown in Table [T] and are the same as in Hirashita & Yan (2009) and Kuo & Hirashita (20ll). 
Although we calculate silicate and carbonaceous dust separately, we are interested in how the overall grain size distribution is affected by 
each dust process. Therefore, we focus on the total grain size distribution. 

In this work, we assume that grains are spherical. Thus, the mass of a grain with radius a is 

m(a) = — g— s, (9) 

where s is the bulk density of dust grains. In our model, we consider that the minimum and maximum radii of grains, amin and Omax, are 
0.0003 [im and 8 /xm, respectively. Although the minimum size of grains is po orly known, eve n if amin = 0.001 /xm, the evolution of both 
the total dust mass and the grain size distribution does not change significantly l lHirashitall20I2h . 



2.2. 1 Dust production by AGB stars 



The size distribution of grains produced by AGB stars is not well known. Winters et al. 1119971) suggested that the size distribution is log- 



normal with a peak at ~ 0.1 ^m based on the fitting to observed SEDs. Yasuda & Kozasa 



( j20I2h 



have recently calculated the size distribution 



of SiC produced by C-rich AGB stars by performing dust formation calculation coupled with a hydrodynamical model. They showed that the 
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mass distribution, f[a), is close to log-normal with a peak at 0.2-0.3 ^m, where the grain size distribution /(a) is defined so that /(a)da 
is the number density of dust grains with radii in the range [a, a + da] (The size distribution multiplied by a* means the mass distribution 
per logarithmic grain radius). Hence, both theory and observations suggest that AGB stars preferentially produce large grains (a > 0.1 /xm). 
In this paper, we simply assume that the mass distribution, a'^f{a), of each species produced by AG B stars is l og-normal with a peak at 
0.1 ^m with standard deviation a = 0.47, so that the shape of the mass distribution in Fig. 7 in .Yasuda & Kozasal 20131 is reproduced. We 
normalize /(a) by 

md{m) = / — a sf{a)da, (10) 
Jo 3 

where md(m) is the dust mass released by a star with mass m. The size distributions of all species are assumed to be the same for sim- 
plicity. Dust mass data for AGB s tars with mass 1-7 Mq and metalhcity Z = (5.0 x 10~^, 0.1, 0.2, 0.4, 0.75, 1.0) Zq is taken from 



Zhukovska. Gail. & Trielofn ( 120080 . The size distributions of dust species other than carbonaceous dust are summed to compose the grain 



size distribution of silicate (the same procedure is also applied in Section |2.2.2| l. We define /x(ffl) as the size distribution of dust species, 
where X represents the key element of dust species (X = C for carbonaceous dust and X = Si for silicate dust). 



2.2.2 Dust production by SNe 11 

Some fraction of dust grains in galaxies are produced in the ejecta of SNe II (e.g.. lMatsuura et al Hoi ij). After a SN explosion, reverse shock 
occurs because of interactions between the ISM surrou nding the SN and its ejecta, and dust grains are destroyed by sputtering in the shock 
(e.g., iBianchi & Schneider 2007. ; .Nozawa et al..,2007h . Nozawa et al. 1 2007 ) calculated the total mass and size distribution of dust grains 
ejected by SNe II considering t he dust destruction in the radiative and non-radiative phases of SN remnants. We adopt the data for dust mass 
and size distribution derived by Nozawa et al. 720 071) for SNe II with mass 13-30 MoB They also considered two cases for mixing in the 
helium core: unmixed and mixed models. Hi rashita et al.l ( 120051) showed that the data from the unmixe d model is in better agreement with 
observations than that of the mixed model. Hence, we adopt the unmixed model. iNozawa et alj ( 120071) showed that the size distribution of 
grains supplied by SNe II is biased to large (~ 0.1 fim) grains due to the destruction of small grains by the reverse shock. 

The amount and size distribution of grains injected by SNe II depend on the density of the surrounding ISM because the dust destruction 
efficiency of the reverse shock is higher in the denser ISM. However, the trend that smaller grains are more easily destroyed does not change. 
In this paper, the hydrogen number density of the ISM surrounding the SNe II, nsN, is set to be 1.0 cm~'^ as a fiducial value, but the cases 
with nsN = 0.1 and 10.0 cm^'^ are also examined. 



2.2.3 Dust destruction by SN shocks in the ISM 



Dust grains in the ISM are destroyed or become smaller by sputtering due to the passage of interstellar shocks driven by SNe. Since the 
destruction changes and depends on the grain size, it is important t o consider dust destruct ion taking into account the grain size distribution. 

To calculate this destruction process, we adopt the formulae in lYamasawa et al.l j201lh . which we show here briefly. The number density 
of dust grains with radii in the range [a,a + da] after the passage of a SN shock, /x(a)da, is given by 



/x(a)da : 



5x(t[, a')dafx{a')da' , 



(11) 



where (,x{a, a')da is the number fraction of grains that are eroded from the initial radii [a' , a' + da'] to radii [a, a + da] by sputtering in 
the SN shock and has been obtained using the models bv lNozawa et alj (|2006}. Note that if a > a', ^x(a, a') = 0. Thus, the change in the 



number density of grains with radii [a, a + da], dNd,x{a-), after the passage of a single SN shock is expressed as 



diVd.x(a) 



^yi{a, a')dafx(a')da' — /x(a)da. 



(12) 



Accordingly, the change of mass density, dMd,x(ci), is 
dMd,x(a) = ^na^sdNdMa) 



4:Tva s 



^x(ffl, a')da/x(a')da' — A^d,x(a)da, 



(13) 



where A^d,x(ci)da = |7ra^s/x(a)da — Md,x{a) is the total dust mass with radii [a, a + da] before the dust destruction. The dust 



^ Although iNozawa et alj )2007ll) investigated only the dust formation in SNe II evolving from zero-meta llicity stars, th e grain species for med in the ejecta 
of SNe II and their size distribution are insensitive to the metallicity of progenitor stars (e.g., IXodini & Ferrara 2001 ; K ozasa et alj|2009h . In addition, the 
destructio n efficiiency of dust by the reverse shocks is almost independent of metallicity in the ISM; its difference between Z = and Z = Zq is less than 
15% (see lNozawa et alj2007l) . 
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destruction efficiency depends on the hydrogen number density of the ISM, nsN, such that dust grains are destroyed more efficiently in 
denser regions. As mentioned in Section l2.2.2l nsN ~ 1-0 cm^'^ as a fiducial value in this paper. 

The equation for the time evolution of Md,x(a, t) for dust destruction by SN shocks in the ISM is expressed as 

dMd,x(a,i) _ Mewcpt 



dt 



MlSM{t) 



7SN(t) 



Md.x{a,t) ^ rn{a) I ^x(ci, a )da/x(a , t)da' 



(14) 



where ^sN{t) is the SN rate and A/swcpt is the ISM mass swept up by a SN shock. To express the dependence on the galaxy age, we write 
Md,x{a) and /x(a) as Md,x(a, t) and /x(a, t), respectively. 
The SN rate, 7sN(t), is expressed as 

(•40 MfT, 



7SN(t) = 



(m)SFR(t-r„)dm, 



(15) 



8 Mq 



where we assume that the range of integration in Eq. (115) is the mass range where SNe occur (IHeger et al .l2003h (ift-r„ < 0,SFR(t 
0). 

The ISM mass swept up by a SN sh ock, Afswopt, depends on the density and metallicity of the ISM. In our model, we adopt the following 



formula used in 

A</swept 



Yamasawa et al. 



1 201 II) 



M, 



1535nq 







+ 0.039 



(16) 



2.2.4 Grain growth 

Here, we formulate the growth process of grains taking into accou nt the grain size distribution. In the ISM, particularly in dense and cold 
regions, metals accrete onto the surface of grains efficiently (e.g.. lLiffman & Clavtonll 19891; llnouell2003l:lDra inell2009i). Recently, vario us 



studies have shown the importance of grain growth for dust enrichment i n galaxies (e.g., Zhukovska. Gail. & Trie loff 2008^; 



Michalowski et al 



ill „ . - ■■ - . 

Valianteet al.ll201ll : IHirashita & Kuol boi lb . Iffirashita & Kuol hoi l|) sho wed quantitatively that the grain size 
distribution has a very important consequence for grain growth. Here, we follow the formulation by Hirashita & Kuol ( 2011 ) and consider 
only grain growth of refractory dust (silicate and carbonaceous dust in this paper). Although volatile grains such as water ice also exists in 
clouds in reality, they evaporate quickly when the clouds disappear or the gas temperature rises. 
For grain growth, the following equation holds: 



9fx{a,t) , d _ 

where d = da/dt is the growth rate of the grain radius 0. 
From Eq. (|9), 

dm(a) 



da 



= 47ra s. 



Also, from Hirashita & KucI 1 2011 ), the rate of mass increase of a grain with radius a is expressed as 
dm(a) 



dt 



(17) 



(18) 



(19) 



where is the mass fraction of the key species X in the grains, mx is the atomic mass of X, a is the sticking coefficient of the key species, 
and IZ is the collision rate of X to a grain with radius a, defined as follows i Evansl[l994 ) 



TZ — A-Ka nx 



1/2 



(20) 



^ 27rmx / 

where n^(t) is the number density of X in the gas phase in the CNM, k is the Boltzmann constant, and Tgas is the gas temperature in the 
CNM. We apply gy^ = 1.0 and 0.166 for carbonaceous and silicate dust, respectively (Table 1) and Tgas = 100 K. 
Next, we consider nx(t), which is evaluated by 



nx(t) = 



mx Mism{t) 



5x- 



Md,x(t) 



Mx(t) 



(21) 



where Pi^M is the average mass density of the ISM in which grain growth occurs. As grain growth occurs, the number of gaseous metals de- 
creases. Thus, nx is a decreasing function of time if only grain growth is concerned. The mass density is estimated as Pj^m = pmnnH.CNM, 
where p is the mean atomic weight, assumed to be 1.33 (the mass ratio of hydrogen to helium is 3 : 1). In addition, mn and nn.CNM are 



the mass of a hydrogen atom and the hydrogen number density in the CNM, respectively, and we apply ncNM = 30 cm 
above four equations [Eq. ( I18M21H . we obtain 



. Hence, from the 



Note that Eq. (TTJ is valid for the case where only grain growth is considered, i.e., without sputtering, shattering and coagulation. 
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da 

dt 



«PlSM 



5xs 



MisM(t) V27rmx/ 



1/2 



5x- 



Mx(t) 



(22) 



We assume a — 1 for simplicity, wiiichi means thiat wiien tiie key species collide with a dust grain, it definitely sticks. In our study, we 
calculate the grain growth using Eqs. ( 117) and ( 122) . 



2.2.5 Shattering 

Turbulence occurs in the I SM ubiquitously, and it is confirmed that turbulence is maintained by thermal conduction from simulations (e.g., 
Kovama & InutsukjEooj) . In a turbulent medium, dust grains are accelerated by turbulence (e.g., McKee & Ostriker 2007 ), and they col- 
lide with each other and shattering can occur (e.g., lYan et aljEooi : Iffirashita & \^l2009l : iHirashita et al .I l2010ll . iHirashita & YanI 1200^ 
suggested that the grain size distributi on in the ISM c hanges significantly by shattering due to collisions between dust grains accelerated 
by magn etohydrodynamic turbulence ( lYan et al .112004ft . In our model, to calcula te shattering pro c ess, w e adopt the grain velocity calcu- 
lated by Yan et al. ( 2004 ), a nd the shattering equation and parameters used by IHirashita & YanI ( l2009t) , whose formulation is based on 



Jones. Tielens & HoUenbaclj jl99d) . 

We outline the treatment of shattering. We define px(a, t)da — m{a) fx{a, t)da as the mass of grains with radii [a, a + da] in a unit 
volume (refer to as "mass density" in this paper). Considering shattering in the collision between two grains with radii ai and a2 (called 
grains 1 and 2, respectively), the time evolution of px{ci, t)da for shattering is expressed as 

dpx(a, t)da 



dt 



m{a)px(a,, t)da I ct[m{a), m{ai)]px(ai,t)daidai 



a[m{ai), m(a2)]px(ai, t)daipx(a2, t)da2mj^^j(a)daida2, 



and 



a[m{ai) , m{a2)] = 





m{ai)m{a2) 



{Vl,2 ^ fahat) 
(«1,2 > ^shat), 



(23) 



(24) 



where 'ml^^^{a) is the total mass of shattered fragments of grain 1 within size bin [a, a + da] by a collision betwe en grains 1 and 2, and de- 



pends on the relative velocity of the grains. The size distribution of shattered fragments is proportional to a (e.g., Jones. Tielens & Hollenbacl] 



ll99 d)FI. The cross section of the collision between grains 1 and 2 is assumed to be cri,2 = 7r(ai + a2)^]- The shattering threshold, t ^ahat. 
is assumed to be 1.2 km s^^, and 2.7 km s^^ for c arbonaceous dust a nd s ilicate dust, respectively ( Jones. Tielens & HoIlenbacl]|l996h . We 
adopt the same treatment for the relative velocity as I Jones et al. ( 1994 ) and Hirashita & Yan ( 2009 ): Each time step is divided into four small 
time steps, and we consider shattering under the following four relative velocities in each small time step (i) front collision (?;i,2 = vi+ V2), 
(ii) back-end collision (111,2 = l^i — 112 1), (iii) side collision ui,2 = vi, and (iv) i;i,2 = V2- Here, iii and V2 are the velocities of the grain 

with radius ai and 02, respectively. 

Shattering can occur not only in turbulence but also in SN shocks (e.g., Jones. Tielens & HoUenbact] 199d) . However, both of these 
shattering mecha nisms have similar consequences on the grain size distribution, so it is difficult to separate them. To compare our work with 
previous studies ( Hirashita et al.ll201ot Kuo & Hirashitall2012 ). we only consider shattering in turbulence. 



2.2.6 Coagulation 

In low temperature and high density regions of the ISM, it is expected that coagulation by grain-grain collisions occurs. Indeed, Stepnik et alj 
1 20030 observed dense filaments and showed that the ratio of the intensity in the filaments, Jeo /jm/ ^100 /jm, is smaller than that in the diffuse 
ISM. They concluded that this trend resulted from the decrease of small gr ains due to coagulation. For coagulation, we adopt the formulation, 
the velocity of grains, and the parameters used bv fairashita & Yan|j2009l) . 
The time evolution of px{a, t)da for coagulation is expressed as follows 



dpx{a, t)da 



dt 



+ 



and 

a[m{ai),m{a2)] = 



m{a)px{a, t)da 



'.[m(ai),m{a)]px{ai, t)daidai 



a[m{ai) , m{a2)]px{ai,t)daipx{a2, t)da2mcis,g{a)daida2, 



/3cri,2«l,2 
m{ai)m[a2) 



ivi,2 ^ vl^lg) 



(25) 



(26) 



^ The method of calculation of the maximum and minimum size of fragments is described in detail in Section 2.3 in iHiiashita & Ya^ i2009h . 
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where /3 is the sticking coefficient of dust grains, and m]:^^^{a) — m{a\) if the mass range of m(ai) + 771(02) is within [m{a),m{a) + 
dm(a)]; otherwise mcoig(a) = 0. 

We assume that coagulation occurs if the relative velocity is less than the coagulation threshold t^ciag- In our model, it is calculated in 
the same way as lHirashita& Yai]j2009h . 



a\ + a\ 
(fli + 



1/2 



(27) 



where _Ri,2 = aia2/(ai + 02), 7 is the surface energy per unit area, and 1/ E — [(1 — vif / Ei + (1 — V2f' / £2], where vi and Ei are 
Poisson's ratio and Young's modulus of grain 1. The parameters we used are shown in Table [T] Here, we assume /3 = 1 for simplicity. The 
treatment of the relative velocity is the same as for shattering. 



2.2. 7 Formulation of the grain-size dependent evolution of dust mass 



Here, using the dust processes introduced above, we show the equation for the dust mass evolution in a galaxy at each grain radius bin, so 
that we can finally obtain the evolution of grain size distribution. Defining AA/d(a, t) = m(a)f[a, t)Aa as the mass density of grains with 
radii [a,a + Aa] □ and a galactic age t, it is formulated as 



dAMi{a,t) 
di 



AMi{a,t) 

MlSM(t) 

-^^swcpt 



+ AYd{a,t) 



7SN(t) 



MiSM(t) 

m{a)Aa- 



+7?CNM 
+'?WNM 
+77WNM 



AMi{a,t) - m{a) 



^{a,a )Aaf{a ,t)da' 



dt 



dAMd{a,t) 



dt 

dAMi{a,t) 



+ riCNM 



shat.WNM 



dAMd{a,t) 



dt 



+ 7?CNM 



coag.WNM 



dt 

dAMd(a,t) 
di 



shat.CNM 



coag.CNM 



(28) 



where 77WNM and 77CNM are the mass fraction of WNM and CNM in the ISM, respectively. From top to bottom, the terms on the right hand 
side describe reduction of dust due to astration, ejection of dust from stellar sources, dust destruction by SN shocks, grain growth in the 
CNM, shattering in the WNM and CNM, and coagulation in the WNM and CNM. To calculate the dust processes which occur in each ISM 
phase, as mentioned in Section [T] we assume (1) that t^wnm and t^cnm are constants and (2) that there are two stable phases, WNM and 
CNM, in the ISM (namely, the sum of rywNM and r;cNM is unity). 

The total mass of grains with radii [a, a + Aa] ejected by stars per unit time, AY'd(a, t), is expressed as 



AYd{a,t) = 



Amd(m, Z{t — r,„), a)(f>{Tn)SFR{t — rm)dm, 



(29) 



"Icut (t) 



where Amd(m, Z, a) is the total mass of grains with radii [a, a + Aa] released by stars with mass m and metallicity Z. 



3 MODEL RESULTS 



In this paper, as mentioned above, we consider the effects of dust formation by SNe II and AGB stars, dust destruction by SN shocks in the 
ISM, grain growth, shattering, and coagulation on the evolution of grain size distribution in galaxies. Among these processes, dust formation 
by SNe II and AGB stars, dust destruction and grain growth directly increase or decrease the total dust mass, while shattering and coagulation 
modify only the grain size distribution. The evolution of the total dust mass in galaxies is often modele d by taking into account the former 
four contributions (dust formation by SNe II and AGB stars, dust dest ruction, and grain gro wth) (e.g., Dw ek & Scalo I980t Dwek I998I : 
Calura. Piping & Matteuccill2008l:Ehukovska. Gail. & Tr ielofl 2008; Pi pino et alj201ll : llnoue 201 1; Asano et alj20I3l) . They calculated the 
dust evolution by assuming a representative grain size, but the dust destruction and grain growth depend on the grain size distribution. Thus, 
it is unknown whether these four contributions can reproduce the grain size distribution in nearby galaxies even though they can explain the 
evolution of the total dust mass. In Section lBTl we first investigate the contributions of the processes that directly affect the total dust mass, 
and then in Section lT2l we examine the effects of shattering and coagulation. 



In this Section, we use the symbol "A" to emphasize that it stands not for infinitesimal but a certain small amount. 
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Figure 1. Grain size distribution talcing into account the dust production by AGB stars and SNe II and dust reduction through astration (solid lines). Red, blue, 
and purple lines represent the cases aX t = 0.1, 1.0, and 10 Gyr, respectively, with rgp = 5 Gyr and TigN = 1.0 cm^'^. Dashed lines are cases with dust 
production by SNe II only and dust reduction through astration, the same color corresponding to the same age. Note that the red dashed line overlaps with the 
red solid line. 




3.1 Without the effects of grain-grain collisions 

3.1.1 Stellar processes 

First, we consider the stellar processes including dust ejection from stars (SNe II and AGB stars) and dust reduction via astration. Figure [T| 
shows the result. The size distribution is expressed by multiplying a'^ to show the mass distribution in logarithmic grain radius bin. We adopt 
rsF — 5 Gyr and tisn ~ 1-0 cm""^. We also show the cases with SNe II only. As mentioned in Section|2] since Aftot is just a scale factor, 
the shape of the size distribution does not change with A/tot- From Fig. [T] throughout any galactic age, we can observe that the grain size 
distribution has a peak at around 0.5 ^m, and that only a small amount of grains with a < 0.01 fim can be formed by stars. In Fig.|2] we 
show the grain size distribution for other values of nsN^ 0.1 cm""^ in the left panel and 10.0 cm""^ in the right panel, respectively. From 
Figure|2] we find that a larger amount of dust grains with radii less than ~ 0.1 ^im are destroyed by reverse shocks in the case of higher nsN, 
and a smaller amount of dust is supplied to the ISM. However, even if nsN changes, the trend that a small amount of dust grains with radii 
less than 0.01 fim are supplied to the ISM does not change. Thus, stars are the sources of dust grains with large radii (^ 0.05 ^m). 
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10 



10^ - 




grain radius [yu,m] 

Figure 3. Grain size distribution taking into account the dust destruction by SN shocks in the ISM in addition to the processes in Fig. [T] (dotted lines). The 
values of rgp and are the same as in Fig.[T] Solid hues are the same as in Fig.[T] the same color coiTesponding to the same age. Note that the red dotted 
line overlaps with the red solid line. 



From Figs. [T] and |2l we observe that dust from SNe II always dominates the grain size distribution, while the contribution of dust from 



AGB stars is seen only around 0.1 /^m at a galactic age t — 10 Gyr. From our calculation, for the case with tisn 

-3 



1.0 cm the dust 



mass ratios produced by AGB stars and SNe II are 1.6 x 10 , 0.16, and 0.37 at t = 0.1, 1.0 and 10 Gyr, respectively. On the other hand, 
Valiante et al] j2009h suggested that the contribution of AGB stars to the total dust mass in gal axies approaches or e xceeds that of SNe II at 



Nozawa et al. 



2007h, whereas .V aliante et al 



t ^ 1 Gyr. This differe nce mainly results from the d ust mass data ado pted. We adopt the da ta of|l _^ 

1 20091) adopted those of lBianchi & Schneideil J2007h . The dust mass of lNozawa et alj j2007h is larger than that of lBianchi & Schneideil ( 12007 
because of the difference in the treatment of the dust condensation and the destruction by reverse shocks. However, even if the contribution 
of AGB stars is larger, the result that only a small amount of grains with < 0.01 /im are produced by stars does not change. 

From the right panel in Fig. |2] we find that the contribution of dust from AGB stars is relatively large for nsN ~ 10.0 cm^'^ at 
t = IQ Gyr. At t = 10 Gyr, the dust mass ratios produced by AGB stars and SNe II are 0.16, 0.37, and 1.39 for the cases with rzsN = 0.1, 1.0 
and 10.0 cm~"^, respectively. This is because a larger amount of dust grains condensed in the ejecta of SNe II are destroyed by reverse shocks 
for higher nsN- 



3.1.2 Dust destruction by SN shocks in the ISM 

In Fig.m we show the evolution of the grain size distribution taking into account dust destruction by SN shocks in the ISM in addition to the 
dust production by SNe II and AGB stars. We also present the cases without the dust destruction (i.e., the same as the solid lines of Fig.[T}. 
The values of tsf and tisn are set to the same values as in Fig.[T] Alt< 1.0 Gyr, the grain size distributions with and without the dust 
destruction by SN shocks in the ISM are very similar to each other. 

Now we estimate the dust destruction timescale. First, we introduce the sweeping timescale, Tswecp, at which SN shocks sweep the 
whole ISM, as 

_ MiSM 
"^swccp — , ^ • 
-/Wswcpt7SN 

From Eq. dlSt . if C is defined as 

/40 Mq 
(^(m)dm, (31) 
. Mq 

Eq. ( list , with Eq. (0, can be approximated as 

7SN ~ C , (32) 

TSF 

where C is about 1.5 x 10~^ from our calculation. Thus, if nsN ~ 1-0 cm~^, Tswecp ~ 2-4 x 10~^tsf- Next, we approximate the dust 
destruction rate by introducing the dust destruction timescale, tsn, as 

dMd ^ _ Md 

at SN TSN 
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grain radius [/um] grain radius [/U,m] 

Figure 4. Same as in Fig. [3] but we adopt different values of nsN^ 0.1 cm~^ in tlie left panel and 10.0 cm~'^ in the right panel. Note that the red dotted lines 
overlap with the red solid lines. 



The right hand side of Eq. \\A\ can be approximated as —r^JccpMdi)^ — £,) where ^ is a typical value of 5x(ci, a ); then, Eq. ( I14t reduces to 

rSN ~ (1 — (,)~^ Ts-^cap ■ (34) 

Since the overall efficiency of dust destruction, (1 — ^), is ~ 0.3 for nsN ~ 1-0 cm~^ I Nozawa et al J200^ . we obtain tsn ~ O.lrsp. Thus, 
the difference between the cases with and without dust destruction cannot be seen at f = 0.1 Gyr in Fig.|3] where tsn ~ O.Itsf ~ 0.5 Gyr. 

We find that dust grains with a < 0.1 /im are destroyed effectively at 10 Gyr ^ tsn (compare the solid and dotted lines in Fig.|3j. 
Since the decreasing rate of grain radius by sputteri ng does not depends on the grain radius, smaller grains are effectively destroyed in SN 
shocks, and the amount of smaller grains decreases ( Nozawa et alj20()^ . 

In Fig.H we show the cases with nsN = 0.1 and 10 cm . Comparing figs.[3]and|4] we find that a larger amount of dust grains 
are destroyed for higher nsN- The destruction effect is more pronounced at small sizes. Indeed, we observe that grains with a < 1.0 /im 
are effectively destroyed in the case with nsn = 10.0 cm~^. Nevertheless, the result that smaller grains are effectively destroyed does not 
change, and we find that dust grains with radii larger than 0.1 /im mainly dominate the total dust amount in galaxies. Consequently, if the 
dust destruction by sputtering in SN shocks is dominant, only large (a > 0.1 fxm) grains can survive in the ISM. 



3.1.3 Grain growth 

Figure |5] shows the evolution of the grain size distribution taking into account the dust production from stellar sources, dust destruction, and 
grain growth. We adopt rsp = 5 Gyr, nsN = 1-0 cm~^, and the mass fraction for the CNM, r/cNM = 0.5. From Fig.|5] we observe that 
while the grain size distributions with and without grain growth are almost the same at ages 0.1 Gyr and 1 Gyr, the difference is clear at 
10 Gyr. The effect of grain growth is prominent around a ~ 0.3 /j,m at 10 Gyr, since the total surface area of grains is dominated by grains 

with a ~ 0.3 /im. The timescale of grain growth is discussed in detail in Section|4l 

In Fig.O we also plot the slope of the MRN size distribution, /(a)da oc a^'^'^da i Mathis. Rumpl & Nordsiec5ll977 ). which is thought 
to be the grain size distribution in the Milky Way. From Fig.|5] it is clear that the small grains with a < 0.01 /im are too few to reproduce the 



MRN size distributi on. However, th e existence of the 70 /im excess is considered to be a proof of the existence of small grains dBemard et al. 



200 sh . Furthermore, Takeuchi et al . (2003, 2005) argued by using their infrared SED model that small grains are necessary to reproduce the 



near-mid infrared SEDs of star forming galaxies. Consequently, when we consider the case in which dust production by SNe II and AGB 
stars, dust destruction, and grain growth take place, the grain size distribution is always dominated by large grains, and we need to consider 
other processes to produce small grains efficiently. 



3.2 Grain-grain collision effects 

In the above we have investigated the dust processes which directly affect the evolution of the total dust mass in galaxies: dust production by 
AGB stars and SNe II, dust destruction by SN shocks, and grain growth. As shown above, these processes cannot produce small grains (a < 
0.01 /im) efficiently. Therefore, we now consider the contributions of the grain-grain collisions, shattering and coagulation in turbulence, to 
the grain size distribution. If these processes occur, although the total dust mass in galaxies does not change, the grain size distribution does. 
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Figure 5. Grain size distribution taking into account the dust production from stellar sources, dust destruction, and grain growth with rsp = 5 Gyr and 
TigN = 1.0 cm~'^ (long-dashed lines). We adopt ??CNM = 0-5- Dotted lines represent cases without grain growth [the same as in Fig. (3j]> the same color 
corresponding to the same age. Note that the red dotted lin e overlaps with the red long-dashed line. Green solid line represents the slope of the power-law grain 
size distribution with index —3.5 [/(a)da cx: a~^'*da iMathis. Rumpl & Nordsiec3ll977h l which is thought to be the grain size distribution in the Milky 
Way. 



(D 
> 



0) 
^1 



CO 
CO 



10' 



10' 



« 10"^ 



10' 



10" 



10 



-4 



\ 

\ 



/}■/ 

.A 
'■'li 
n 

■ M . . 



1/ ^ 

fl^ Red : O.lGyr \| 
Blue : IGyr \ V 



Purple: lOGyr 



10"'* 10"^ 10"^ 10"' 



10° 10' 



grain radius [/xm] 



10.000 

£i 1.000 



o 
o 



s 

■a 



0.100 



0.010 - 



0.001 



w/ shattering 
w/o shattering 

w/o grain growth 
& shattering 



metallicity 



0.1 1.0 10.0 

Galactic age [Gyr] 



Figure 6. Left panel: grain size distribution with (dot-dashed lines) and without (long-dashed lines) shattering (all other dust processes in Section [Til are 
included). Note that the red dot-dashed line overlaps with the red long-dashed Une. Right panel: time evolution of dust-to-gas mass ratio with (solid line) 
and without (dotted Une) shattering. The case without grain growth and shattering (dashed line) is also plotted. Dot-dashed line represents the evolution of 
metallicity. The parameters rsp and nsN are set to be 5 Gyr and 1.0 cm^^, respectively. We adopt rywNM = ^?CNM = 0.5. 



3.2.1 Shattering 

In the left panel of Fig. |6l we show the evolution of the grain size distribution in the galaxy with and without shattering (all other dust 
processes in Section lSTI are included). The right panel of Fig.|6]shows the time evolution of dust-to-gas mass ratio (Afd/AfisM) for the cases 
with and without shattering, respectively. We also plot the case without grain growth and shattering and the evolution of metallicity in the 
same panel. The parameters tsp and nsN are set to be 5 Gyr and 1.0 cm"'', respectively. We adopt j7wnm ~ Vcnm = 0.5. 

From the left panel of Fig.|6] at the early stage of galaxy evolution (0.1 Gyr) the size distributions with and without shattering are 
similar with only a little difference at small sizes. At 1 Gyr, we observe that the size distribution has a bump at a ~ 0.001 ^im in the case 
with shattering. As time passes, the amount of large grains decreases, and as a result the size distribution is dominated by small grains. This 
behavior is substantially different from that of the case without shattering. We now discuss this behavior in more detail. 
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Figure 7. Left panel: grain size distribution with (triple-dot-dashed lines) and without (dot-dashed lines) coagulation (all the other dust processes are included). 
Note that the red and blue dot-dashed lines overlap with the red and blue triple-dot-dashed hues. Right panel: time evolution of dust-to-gas mass ratio with 
(solid Une) and without (dotted line) coagulation. Dot-dashed line represents the evolution of metallicity. The values of parameters (rgp, ngt-j, JTwNMi and 
?7CNm) are the same as in Fig.|6] 



As shown in Eq. ( 1231 ), the efficiency of shattering is larger for larger amount of grains (IHirashitalEOlOb . At 0.1 Gyr, the efficiency 
of shattering is low because of the small dust abundance. As a result, there is only a small difference between the cases with and without 
shattering. Alt = 1 Gyr, since shattering occurs efficiently due to the increased amount of large grains, the amount of small grains increases. 
At the same time, we observe that the grain size distribution has little difference between the cases with and without shattering at a > 
0.1 iim in the left pan el of Fig. [S] This is because shattering of a tiny fraction of large grains can produce a large amount of small grains 
i HirashitaetaP boio'). I^urthermore, since the number of small grains increases, the small grains dominate the total grain surface area. 



Consequently, grain growth occurs at the smallest grain sizes (a < 10"'^ /^m), forming a bump at 10~'^-10~^ ^m. At t = 10 Gyr, since 
the number of small grains increases, large grains are shattered more efficiently by the frequent collisions with the small grains. Consequently, 
comparing the grain size distribution at 10 Gyr with that at 1 Gyr, the amount of large grains decreases significantly. Furthermore, because 
of grain growth, the bump is shifted to a larger size at 10 Gyr than at 1 Gyr, and finally the size distribution has a large bump at a ~ 0.01 /xm 
at 10 Gyr. 

Focusing on the grain size distribution at 10 Gyr, we find that if shattering occurs, the amount of grains with a > 0.2 fj,m is more than 
two orders of magnitude smaller than that of grains with a < 0.2 ^m. Thus, the maximum size of grains in diffuse ISM is determined not by 
Stardust but by the process of shattering. 

In the right panel of Fig. |6] we find that grain growth starts to increase the total dust mass at around t = 1 Gyr as seen in the rapid 
increase of dust-to-gas mass ratio, and grain growth becomes more rapid in the case w ith shattering than in the case without shattering 
because of the increased number of small grains. As discussed in IKuo & Hirashital l2012h , shattering contributes not only to the evolution of 
the grain size distribution but also to the total dust mass in galaxies indirectly through the enhanced grain growth. Thus, shattering is a very 
important process in understanding the evolution of the size distribution and the amount of dust grains in the ISM. 



3.2.2 Coagulation 

In Fig.[71 we show the evolution of the grain size distribution with and without coagulation (all the other dust processes are included) in the 
left panel, and the time evolutions of dust-to-gas mass ratio with and without coagulation and of metallicity in the right panel. The parameters 
adopted are the same as in Fig.|6] 

From the left panel of Fig. |7] we find that there is little difference between the cases with and without coagulation at 0.1 and 1.0 Gyr. 
Since larger grains are coupled with the larger-scale turbulence, they can obtain larger velocity dispersions. Thus, coagulation mainly occurs 
by collisions between small grains whose velocity dispersions are smaller than the coagulation threshold [Eq. (I27H . However, since the 
abundance of small grains is low, the contribution of coagulation is not seen at 0.1 and 1 Gyr before shattering becomes effective. After that, 
a large abundance of small grains are produced by shattering so coagulation becomes effective. Consequently, the bump at a ~ 0.01 /im 
shifts to a larger size by coagulation. 

From the right p anel of Fig . IT] we find that the evolution of the total dust mass does not change significantly by coagulation, confirming 



the result obtained by iHirashital ( 120 1 21) . If coagulation occurs, the number of small grains decreases; as a result, the surface-to- volume ratio 
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grain radius [/U-m] grain radius [/U,m] 

Figure 8. Grain size distribution with (rjwNMi '?CNm) = (0-9: 0-1) (solid lines with filled circles in the left panel) and (0.1, 0.9) (solid lines with filled 
diamonds in the right panel). The values of parameters rsp and nsN are the same as in Fig.|7] Triple-dot-dashed lines are the case with »)wNM = ^CNM = 
0.5. 

of grains decreases. This effect may suppress the increase in dust mass due to grain growth. However, since grain growth becomes inefficient 
to the dust evolution before coagulation becomes efficient (the details are shown in Section |4l, the contribution of coagulation cannot be 
observed for the total dust mass evolution. Consequently, the effect of coagulation on the evolution of the total dust mass in galaxies is 
negligible. 

From the left panel of Fig. [7] we find that the amount of grains with a > 0.2 /^m does not change significantly by coagulation because 
coagulation cannot occur by coUision between large grains which have larger velocity dispersions than the coagulation threshold. Thus, 
although the bump of the grain size distribution is shifted to a larger size by coagulation, coagulation does not affect the maximum size 
determined by shattering (Section [3.2.1t . 

3.3 Parameter dependence 

Shattering and coagulation occur differently in both ISM phases (WNM and CNM). Here, by adopting (77WNM, t;cnm) = (0.9, 0.1) and 
(0.1, 0.9), we show the effect of ISM phases. 

In Fig.[8l we show the evolution of the grain size distribution with (?;wnm, ?7cnm) = (0.9, 0.1) (left panel) and (0.1, 0.9) (right panel). 
The case with r;wNM = ?7cnm = 0.5 is shown for comparison. At t = 0.1 Gyr, the grain size distributions are almost the same in all the 
cases, since the dust process is dominated by the production by stellar sources. At t = 1 Gyr, the amount of dust grains with a < 0.01 /im 
is larger for the case with larger rycNM because grain growth is more efficient. At 10 Gyr, the difference is clear at each grain size. For 
larger rywNM, the decrements of the amount of dust grains with a > 0.2 /im is larger because shattering in the WNM is more efficient. 
Furthermore, for larger ?7cnm, the bump produced by grain growth around 0.01 fim shifts to a larger size. Thus, we understand that the 
amount of dust grains with a > 0.2 jim and the shift of the bump around 0.01 /im are dominated by shattering in WNM and coagulation in 
CNM, respectively. In addition, comparing the two panels in Fig. 8, we find that the dust amount at a '--^ 0.1-0.2 jim tends to be smaller for 
a larger t/cnmEI It means that the amount of dust grains with a ~ 0.1 /im is dominated not by shattering in WNM but shattering in CNM. 
Hence, the grain size distribution in galaxies is finally dominated by processes in WNM for large grains (> 0.2 /im) and by processes in 
CNM for small grains (~ 0.1 /im). 



4 DISCUSSION 

In Section[3] we showed the evolution of the grain size distribution in galaxies for a variety of mixture of dust processes. We found that the 
grain size distribution is dominated by large grains produced by stars (SNe II and AGB stars) in the early stage of galaxy evolution, but as 
the time passes the number of small grains increases due to shattering, and the small grains grow to larger grains by grain growth. After that, 
the size distribution shifts to larger sizes due to coagulation. Thus, we conclude that, while the grain size distribution in galaxies is controlled 

^ The mass ratio of grains with a ~ 0.1-0.2 /xm for the cases between ('»/wNMi '?CNm) = (0-5i 0-5) ™d (0.9, 0.1) is about 0.6 at t = 10 Gyr. 
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by stellar processes in the early stage of galaxy evolution, the main driver to change the size distribution is replaced with the processes in 
the ISM (shattering, coagulation, and grain growth) at the later stage of galaxy evolution. These processes (shattering, coagulation, and grain 
growth) have timescales dependent on the grain size distribution. In this Section, by adopting representative grain radii, 0.001 /im, 0.01 /im, 
0.1 /im, and 1.0 /im, we discuss the evolution of the grain size distribution more quantitatively. 

Figure |9] shows the time evolution of the ratio between the size distribution functions, /(a), and /(a)star, the latter being obtained by 
considering only the stellar processes [the first and second terms in the right hand size of Eq. (|28H. Panels (a), (b), and (c) are the cases with 
(»7wNM, »7cnm) = (0.5, 0.5), (0.9, 0.1), and (0.1, 0.9), respectively, for tsf = 5 Gyr and nsN = 1.0 cm"''. From Fig. |9l we find that 
behavior of /(a)//(a)star depends strongly on the grain radius. As mentioned above, this is because each dust process works at different 
grain radii on different timescales. First, we discuss the evolution of the grain size in panel (a). We find that /(a) / /(a)star at a = 0.001 /im 
starts to deviate from unity at the earliest galactic age among all four grain sizes. The process that causes this increase is shattering. These 
small grains are produced by shattering between large grains produced by stars. Furthermore, as mentioned in Section |3.2.1| since shattering 
of a small number of large grains can produce a large number of small grains, we cannot see the change of /(a)//(a)star for a = 0.1 and 
1.0 /im. 

At f ~ a few hundreds of Myr, /(a)//(a)star at a = 0.01 /im increases. This increase is also due to shattering. The reason why the 
effect of shattering appears at a = 0.01 /im later than at a = 0.001 /im is that the size distribution of shattered fragments is proportional 
to a~^'^ (see Section [2.2.5b . In other words, the shattered fragments become dominant at smaller sizes on shorter timescales than at larger 
sizes. 

At t ~ 1 Gyr, we find that the increase of /(a)//(a)star at a = 0.01 /im is accelerated. This indicates that another process becomes 
efficient, and it is grain growth. As seen from Fig. |6] as grain growth becomes efficient around 1 Gyr, the amount of grains with less than 
a ~ 0.01 /im increases significantly. 

At t ~ 2 Gyr, we find that /(a) //(a)star decrease at all sizes. These decreases are due to coagulation for small grains (a = 0.001 and 
0.01 /im) and shattering for large grains (a = 0.1 and 1.0 /im). As we showed in Section [3.2.21 coagulation mainly occurs between small 
grains. Thus, the coagulation effect cannot be seen at early phase of galaxy evolution when the abundance of small (a < 0.01 /im) grains 
is small. Shattering can also occur effectively if there is a large amount of small grains because of a high grain-grain collision rate with 
small grains [cf. Eq. (I24H . In addition, the main reason why the decrements of grains with a = 0.1 and 1.0 /im are different is shattering 
in different ISM phases. As shown in Section [53] grains with a > 0.2 /im are mainly dominated by shattering in WNM, while grains with 
a ~ 0.1 /im are dominated by shattering in CNM. In summary, at early phase of galaxy evolution {t < 10 Myr), the size distribution is 
dominated by dust grains produced by stars, after t > 100 Myr, the dust processes in the ISM begin to affect the size distribution at small 
size, and at t ~ 2 Gyr (for tsf = 5 Gyr), various dust processes in the ISM affect all sizes of grains. 

Panels (b) and (c) in Fig.|9]show the cases with (77WNM, 77cnm) = (0.9, 0.1) and (0.1, 0.9), respectively. Compared with panel (a), we 
find that /(a)//(a)star at a = 0.01 /im does not decrease at 10 Gyr in panel (b). This is because the timescale of coagulation becomes 
longer for smaller r/cNM. From panel (c) (ttwnm ~ 0.1, r/cNM = 0.9), we find that the decrement at a = 1.0 /im is smaller than those in 
the cases of panel (a) and (b). This is because the efficiency of shattering in WNM is smaller for smaller rywNM. However, from all the three 
panels, we can observe that the timing at which /(a)//(a)star at all sizes changes due to the dust processes in the ISM (in this case, it is 
about 2 Gyr) does not vary significantly by the change of (ttwnm, ?7cnm) for the same star formation timescale. 

In order to discuss the effect of tsf on the size distribution, the results are shown for the same values of the parameters as in the panel 
(a) of Fig. [21 but for tsf ~ 0.5 Gyr in panel (d) and tsf ~ 50 Gyr in panel (e). Compared with panel (a), we find that /(a)//(a)star 
change at earlier stages for shorter tsf at all sizes. This is explained as follows. If tsf is short, the amounts of dust and metals released 
by stars are large at early phases of galaxy evolution. The times cales of shattering and coagulation are inversely proportional to the dust-to- 



gas mass rati o (e.g., Hi rashita 201C; Hirashita & Omukai 2009^, and the timescale of grain growth is inversely proportional to metallicity 



(e.g., lAsano et al. 2013). Thus, for shorter tsf, dust processes in the ISM (grain growth, shattering, and coagulation) begin to affect the size 
distribution at earlier stages of galaxy evolution (~ 0.6, 2 and 5 Gyr for tsf = 0.5, 5, and 50 Gyr, respectively). The timescale of the change 
of ./ (a)//(a)star is roughly estimated to be ~ 1 (rsF/Gyr)^/'^ Gyr (Appendix IbI. We conclude that the grain size distribution in galaxies 
changes drastically through the galaxy evolution because different dust processes operate on the grain size distribution at different ages. 



5 CONCLUSIONS 

We constructed a dust evolution model taking into account the grain size distribution in a galaxy, and investigated what kind of dust processes 
dominate the grain size distribution at each stage of galaxy evolution. In this paper, we considered dust formation by SNe II and AGB stars, 
dust destruction by SN shocks in the ISM, grain growth in the CNM, and grain-grain collisions (shattering and coagulation) in the WNM 
and CNM. 

We found that the grain size distribution in galaxies is dominated by large grains produced by stars in the early stage of galaxy evolution, 
but as time passes the size distribution is controlled by processes in the ISM (grain growth, shattering, and coagulation) and the age at which 
these ISM processes enter depends on the star formation timescale, as ~ l(rgF/Gyr)^''^ Gyr. While dust production by SNe II and AGB 
stars, dust destruction by SN shocks, and grain growth in the CNM directly affect the total dust mass evolution, we found that the grains 
are predominantly large (a ~ 0.2-0.5 /im) and only a small amount of small grains (a < 0.01 /im) are produced by these processes. If 
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Figure 9. Time evolution of the ratio between tlie size distribution functions, /(a), and /(a)atar, tlie latter being obtained by considering only stellar 
processes [the first and second terms in right hand size of Eq. )28H . Panel (a), (b), and (c) are the cases with (jjwNM , ^CNm) = (0.5,0.5), (0.9,0.1), 
and (0.1, 0.9), respectively with rgp = 5 Gyr and ngN = 1.0 cm~^. Panel (d) and (e) are the case with rgp = 0.5 Gyr and 50 Gyr, respectively with 
('TWNM, '?CNm) = (0.5, 0.5). Solid, dotted, dashed, and dot-dashed lines represent the ratio of a = 0.001 /^m, 0.01 fim, 0.1 /jm, and 1.0 /im, respectively. 
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we take shattering and coagulation into account, the grain size distribution is modified significantly by these two processes. In particular, 
shattering indirectly contribute to the large increase of the total dust mass: After small grains (a < 0.01 /im) are produced by shattering, 
grain growth becomes more effective because of the enhanced surface-to-volume ratio. Furthermore, grain growth produces a large bump 
in the grain size distribution around a = 0.01 /im. The effects of shattering in WNM and CNM on the size distribution appear at different 
grain radii: While grains with a > 0.2 /^m are mainly shattered in WNM, shattering in CNM affects grains with a ~ 0.1 /^m. Furthermore, 
the effect of shattering, in particular shattering in WNM, is large enough to determine the maximum size of grains in the ISM. Coagulation 
occurs effectively after the abundance of small grains is enhanced by shattering, and the grain size distribution is deformed to have a bump 
at a larger size (a ^ 0.03-0.05 /^m at t ~ 10 Gyr) by coagulation. We conclude that the evolution of both the total dust mass and the grain 
size distribution in galaxies are related strongly to each other and the grain size distribution changes drastically through the galaxy evolution. 
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galactic age [Myr] grain radius [/U,m] 

Figure Al. Left panel: the star fonnation histoiy with n = 1.0 (sohd Une) and 1.5 (dotted line). Right panel: grain size distribution with n = 1.0 (triple-dot- 
dashed line) and 1.5 (solid line with triangles). We adopted rgp = 5 Gyr, ng^ = 1.0 cm~^, and f?wNM = ^CNM = 0-5 in these plots. Note that the red 
and blue triple-dot-dashed lines overlap with the red and blue solid lines with triangles. 



APPENDIX A: EXAMINATION OF PARAMETER DEPENDENCE 

In this Appendix, we show dust evolution models with parameters different from the values adopted in the main text. 



Al The Schmidt law index n — l.b 

In Fig. lAll we show star formation history (SFH) and the evolution of the grain size distribution with the Schmidt law index n — l.Q and 
1.5. To compute the SFH and grain size distribution by using star formation rate with the Schmidt index n = 1.5, the SFR with the Schmidt 
law index n = 1.5 (SFR1.5) is expressed as 

SFR1.5W = MM^^ (Al) 

where 1/1.5 is a constant. We define the value of i^i.s so that it satisfies the following equation at t = 0: 

Mismjt) 

TT — T^SF- (AZ) 

SFR1.5W 

This is set to compare it with rsp forn = 1 easily. Thus, we obtain 

vi.^ = rsFM°i. (A3) 
From Fig. lAll we find that the results are not significantly different between the cases with the n = 1.0 and 1.5. 



A2 The index of the Salpeter IMF (7=1.35 

Figure lA2l shows the SFH and the evolution of the grain size distribution with the power-law index of the Salpeter IMF q = 1.35 and 2.35 
(fiducial value in this paper). We observe that SFRs are almost the same, but the grain size distributions are different. If q is small, that is, 
a large number of SNe II are produced, the abundance of dust and metals increase earlier than the case with large q. As a result, the dust 
amount of each size of grains [the values of a*/(a)] with q — 1.35 is larger than the case with q — 2.35. The dust processes in the ISM also 
become effective earlier because of the larger dust abundance. However, we find that the trend of the evolution of the grain size distribution 
(at early phases, stars are dominant sources of dust, as time passes, the processes in the ISM become important) does not change. 



APPENDIX B: TIMESCALE OF THE CHANGE OF F{A) /F(^)star 

In Section [4] we found that the timescale of the change of ,f (a)//(a)star of all sizes of grains depends on star formation timescale, and 
the change are due to coagulation for small grains and shattering for large grains. Since both of shattering and coagulation are coUisional 
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Figure A2. Left panel: the star fonnation histoiy with q = 2.35 (solid line) and 1.35 (dotted line). Right panel: grain size distribution with q = 2.35 
(triple-dot-dashed lines) and 1.35 (solid lines with filled squares). We adopted rgp = 5 Gyr, rigN = 1.0 cm~^ and '»/wNM = 'ycNM = 0-5 in these plots. 



processes, the timescales scale with the grain abundance in the same way. In order to evaluate the dependence on the star formation timescale, 
we compare the contributions of stars and shattering. 

First, we consider the stellar contribution [Eq. ( 1291 1. If D is defined as 

(■oo i-lOO M0 

D= / Amd{m, Z{t ~ Tm), a)4>{m)dmda, (BI) 

"'0 J mcut (t) 

with Eq. the stellar contribution can be approximated as 



dMd 



dt 



~ D- 

star TSF 



Then, we consider the timescale of shattering, Tshat- Since shattering is a collisional process, Tghat can be represented as 
1 



(B2) 



(B3) 



'L'TT-grains 

where (a^) is the 2nd moment of a grain size a, v is the relative velocity of grains, and n^rain is the number density of grains, which is given 
by 



4 , 3x Mi 
-7r(a )sngrain ~ f^nu,BhiitmnY7 > 

o AllSM 



(B4) 



where (a'^) is the 3rd moment of a grain size, s is the bulk density of dust grains, nn.shat is the hydrogen number density in the region 
where shattering occurs, and mn is the mass of the hydrogen atom. We assume the contribution of shattering to the amount of dust grains as 
Afd/rshat, and comparing this equation with Eq. ( IB2b . we obtain the relation between shattering timescale and star formation timescale, 

Md 1 



Tshat — TSF- 



MisM D ■ 

In addition, by substituting Eqs. JB3l l and JB4b into Eq. JB5b . we obtain 



"^shat 



|7r(a3)s 



1/2 



(B5) 



(B6) 



To evaluate this value, we adopt s = 3.0 g cm ^, v — 20 km s"^ and KH.shat = 0.3 cm (WNM) as a representative value. Also, from 
our calculation, D ~ 10~^, and {a'^)/{a^) ~ 10~^ cm for dust grains produced by stars. Then, we finally obtain 



r.hat ~ 1 (^gj [Gyr]. 

1/2 

Thus, we conclude that the timescale of shattering, that is, the timescale of the change of f{a) //(a)star, is proportional to rgp . 
This paper has been typeset from a TgX/ KTgX file prepared by the author. 
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